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. The current noise spectrum of a single-electron transistor (SET) coupled to a nanomechanical 

' resonator is calculated in the classical regime. Correlations between the charge on the SET island and 

^ the position of the resonator give rise to a distinctive noise spectrum which can be very different from 

fSJ . that of the uncoupled SET. The current noise spectrum of the coupled system contains peaks at both 

the frequency of the resonator and double the resonator frequency, as well as a strong enhancement 
of the noise at low frequencies. The heights of the peaks are controlled by the strength of the 
, coupling between the SET and the resonator, the damping of the resonator, and the temperature 

' of the system. 

m ; 

I — I. INTRODUCTION 

' Nanoelectromechanical systems in which a nanomechanical resonator is coupled to a mesoscopic conductor, such 
, as a single-electron transistor (SET), form a new class of mesoscopic system in which there is a fascinating interplay 
^ between the electrical and mechanical degrees of freedom. Coupling to nanomechanical degrees of freedom can modify 
Ch , the transport properties of mesoscopic conductors substantially, giving rise to a number of novel phenomena such as 
. ■ electron shuttlingii^ and phonon blockade effectsi^ Furthermore, nanoelectromechanical systems have important 
"j^ ' potential applications as ultra-sensitive force detectors.«S' 

A nanomechanical resonator coated with a thin metal layer can be coupled to a SET as a mechanically compliant 
I voltage-gatei^ Under such circumstances the tunnelling rates of electrons through the SET island, and hence the 
, current flowing, depend sensitively on the position of the resonator iSiSiS It has been suggested that the SET could 
Ch ' act as a quantum-limited displacement detector for micron-sized resonatorsSiiS thus making it potentially very 
O \ important both for read-out of the resonator's motion in force sensing applicationsii and for exploring quantum effects 
I in electromechanical systemsi^Siii However, fluctuations in the charge on the SET island exert a stochastic force back 

\ on the resonator, affecting its motion, and hence also the average current through the SETil2ii'* Furthermore, the 
coupling between the resonator and the SET also has an important influence on the SET current noise. 
^ ' The effect on the dynamics of a nanomechanical resonator of coupling to a mesoscopic conductor has been inves- 
tigated by a number of authorsi2iiiiiSiiS4iLi& For systems where the nanomechanical resonator forms a compliant 
' voltage gate adjacent to either a SETii^iii or a quantum point contactji^ it has been found that the interaction with 
. the electrons drives the resonator into a steady-state very similar to a thermal state, with an effective temperature 
•■ and damping constant. In particular, when the resonator is treated as a quantum system it is found that whenever 
^ the energy associated with the bias voltage applied across the conductor is much greater than the energy quanta of 
' the resonator, the electrons rapidly dephase the resonator and heat it to an effective temperature proportional to the 
bias voltageii^iiSiii An analysis of the dynamics of a resonator coupled to a SET where the resonator is treated as a 
classical system leads to very similar conclusions in the regime where the applied bias voltage is relatively large and 
' the charge dynamics of the SET is adequately described by the orthodox modeh— 

, Although the dynamics of a nanomechanical resonator coupled to a mesoscopic conductor has been investigated for 
(— ( ' a number of different systems and regimes, there has been little systematic study of the concomitant current-noise 
Q \ spectrum of the mesoscopic conductor. Knowledge of the current-noise spectrum in nanoelectromechanical systems 
O is of interest both because it sets limits on the sensitivity with which the current can be used as measure of the 
^ motion of the resonator and, more generally, because it has long been recognized that the current noise of mesoscopic 
• ^ . conductors provides important additional information, beyond that available from simple average-current properties, 
/\ ' about the interactions which the conduction electrons undergo.-? In this paper the current noise of a SET coupled to 
a nanomechanical resonator is calculated in the classical regime using an extension of the master equation approach 
developed in Ref. Q. It is found that the current noise spectrum of the SET-resonator system can be very different 
from that of the uncoupled SET. In particular, the interactions between the resonator and the transport electrons can 
strongly enhance the SET current noise at low frequencies and give rise to additional peaks in the noise spectrum at 
the resonator frequency and at twice that value. 

This paper is organized as follows. In Sec. II the dynamical model of the coupled resonator-SET system that forms 
the basis of the following calculations is introduced. Then in Sec. Ill, the relationship between the overall current 
noise spectrum of the SET and the spectrum of the charge noise on the island and the noise spectra of the tunnel 
currents between the island and the leads is described. In Sees. IV and V the noise spectra of the tunnel currents 
through the SET junctions and the charge noise on the SET island are obtained, respectively. The full current noise 
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FIG. 1: Circuit diagram for the coupled SET-resonator system. The SET island lies between the two tunnel junctions which 
are assumed to have the same capacitance Cj. Under the drain-source voltage shown, Vds, the electrons tend to flow through 
the SET island from right to left. The resonator is situated adjacent to the SET island and acts as a gate capacitor. 



spectrum of the coupled SET-resonator system is obtained in Sec. VI. Finally, Sec. VII contains the conclusions. 



II. MODEL OF THE COUPLED SET-RESONATOR SYSTEM 

A schematic diagram of the SET-resonator system is shown in Fig. ^ For simplicity, we have assumed a symmetric 
arrangement of capacitances and voltages. The SET island is connected to current leads by tunnel junctions with 
capacitance Cj and resistance R. The total voltage applied across the SET island is Vdg. The nanomechanical resonator 
forms the gate capacitor, with a capacitance Cg{x), and a voltage Vg applied to it. The gate capacitance depends 
on the displacement of the resonator, cc, about its equilibrium separation from the SET island, d. The resonator is 
displaced from its equilibrium position by fluctuations in the SET charge and by thermal fluctuations of the resonator 
itself, but for realistic systems^ ''' the resulting motion is on a scale much less than the equilibrium separation (i.e. 
X <C d), hence an expansion of the capacitance to linear order is sufficient and we may write Cg{x) — Cg{l—x/d), where 
Cg is a constant AiSii^ The gate capacitance is also assumed to be much less than that of the junction capacitances, 
Cg <C Cj, which reflects the situation in typical devices. 

Under a relatively wide range of conditions, the coupled dynamics of a nanomechanical resonator and a SET can be 
described using a master equation formalismi^ Within this formalism, the SET is assumed to be within the regime 
where the orthodox model is valid and the resonator is treated as a classical harmonic oscillator, with angular frequency 
uq, and effective mass m. The effects of damping processes on the resonator, beyond those due to interactions with 
the SET, can be taken into account by including an extrinsic damping rate, 7e. The orthodox model describes electron 
transport through the SET away from the Coulomb blockade regions so that processes other than sequential tunnelling 
of electrons may be neglected.^" The nanomechanical resonator is expected to behave classically whenever the bias 
voltage applied to the SET is much greater than the energy quanta of the resonator, i.e. eVds ^ fiuio}^ Complete 
details of the derivation, justification and analysis of this classical model for the coupled SET-resonator system are 
given m Ref.IIl 

When the background temperature of the SET-resonator system, T^, is much less than the charging energy of the 
SET island, e^/2Cj » fceTe, and the bias voltage is not too large, the charge state of the SET island can be assumed 
to be limited to two possible values, N and N + 1 excess electrons, and a pair of coupled master equations can be 
derived for the SET-resonator system. The bias voltage is also assumed to be much greater than the background 
temperature, eVds S> k-^Te, and so the effects of thermal fluctuations on the electron tunnel rates are neglectedii^iSS 
Writing Pfq(^isi+i){x^ % t) the probability of having N{N + 1) electrons on the SET island and finding the resonator 
at position x and with velocity m, the master equations for the evolution of the probability distributions take the 
formii 

dPw+i 2/ .OPn+i dPN+i d ( G dPw+A „ ^ „ thujIxq \ 
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where P{x,u;t) = PN{x,u\t) + PN+i{x,u;t), G = 2mkBTeje and ^l{r) is the position-independent part of the 
tunnel rate through the L[R) junction. The length-scale Xq has a simple interpretation: it is the amount by which 
the equilibrium separation between the SET and the resonator changes when the charge on the SET island changes 
from N to N +1. 

The master equations [Eqs. and ^] assume that the bias voltage applied across the system is the dominant 
energy-scale and the probability distributions for the resonator are strongly peaked at a; a;o. Under such cir- 
cumstances, which are readily met in practicable resonator-SET systemSj^l^ tunnelling processes that go against the 
direction imposed by the bias voltage may be neglected. The resulting master equations can be used to extract 
equations of motion for the average properties of the resonator and SET, and to characterize the steady-state of the 
system. Furthermore, the master equations provide the necessary information about the coupled dynamics of the 
SET-resonator system to allow the derivation of the SET current noise spectrum. 

Analysis of the coupled master equations^"' shows that the SET electrons act on the resonator like another thermal 
bath (in addition to the substrate to which the resonator is attached) which can be characterized by an intrinsic 
temperature, Ti, proportional to Vds, and an intrinsic damping constant, 7^. The overall effective temperature of the 
resonator is then determined by an average of the temperatures of the two baths at temperatures Ti and Tg, with 
weights that depend on the respective damping constants 7^ and 7e-^^ 

The strength of the interaction between the resonator and the SET is conveniently described by the dimensionless 
parameter k — TOWQa;g/(eVds) and the frequency by e = wot^, where Tt = Re/Vds is the electron tunnelling time. 
In terms of these dimensionless parameters, the intrinsic damping rate of the resonator due to interactions with the 
electrons isi^ 7^ — Ke^/rt- The interaction of the resonator with the electrons also leads to a slight renormalization 
of resonator frequency, e' = e(l — k)^^"^. Finally, the external temperature is conveniently expressed in terms of the 
drain-source voltage as the dimensionless parameter, Q = fceTe/eVds. 



III. SET CURRENT NOISE 



The noise in the current through a SET is generated by the fluctuations in the time dependent current flowing 
through the system. Since the SET contains more than one tunnel junction, the time-dependent current which flows 
is not simply the tunnel current across either one of the junctionsii^iSiiSSiS^ When an electron tunnels through cither 
junction, charge also flows in the external circuit connected to the SET island (so that the charge stored in each 
capacitor returns to its equilibrium value) and hence also contributes to the current. The current flowing through the 
SET can be measured in either the left or right-hand lead and the average current flowing through the leads must 
necessarily be the same. 

The time dependent current through the gated SET is, in principle, slightly different in the left and right-hand 
leads, even for the symmetrical case we consider here where the tunnel-j unction capacitances, Cj, are the samei^ The 
time dependent current through the left-hand lead can be writteni2i2i 

/(^)(i)=a(^)QL(t) + 6(^)Qfl,(t), (3) 

whilst the time dependent current through the right-hand lead is given by2i 

/W(i) = a(^)Qi(t) + fe(^)Qfl(i), (4) 
where Ql(r) are the tunnel currents through the L{R) junctions respectively, andSi 

= = 2c-IW- 

Notice, however, that the two coefficients a*^*) and i = L,R, add up to unity, whichever lead we consider. Because 
of the dependence of these coefficients on the gate capacitance, Cg{x), there will necessarily be a weak dependence 
of the coefficients on the position of the resonator. However, even the leading order position dependent term in 
the coefficients will be smaller by a factor Cg/Cj than those arising within the expressions for the time dependent 
currents. Since we are treating the case where Cg/Cj <C 1, the corrections generated by the position dependence of 
a*^'^ and i — L,R, will be neglected herei^ In fact, it turns out that the noise spectrum is relatively insensitive 
to the exact values of the coefficients, and hence which lead the current is measured in. Therefore, for simplicity, we 
will not distinguish between the two leads in what follows and set a'^' = a'^^^ = a, and 5^^^ = 5^^^ = b. 
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The current fluctuations in the SET are constrained by the conservation of charge. Charge conservation imphes 
that the charge that has tunnelled through the left-hand junction in time t, QL{t), and the charge that has tunnelled 
through the right junction in the same period, Qnit) must be connected by the relation 



At)-QR{t)^QD{t) 



(7) 



where Quit) is the charge dwelling on the island between the two barriers at time t. Furthermore, since the average 
tunnel currents through the left- and right-hand junctions are the same, they are both also equal to the average 
current through the SETji^ 



(8) 



and hence we also have, {Qd) = 0. 

The spectrum of current noise in the SET is given by the Fourier transform of the current-current correlation 
function;^ 



where 



/oo 
dTKii{T) cos(a;T), 
-oo 



Ku{T) = {I{t + r)I{t))-{I{t)f 



(9) 



(10) 



The current-current correlation function, Kii{t), is taken to be independent oft. 

The noise in the total current is a combination of the noise in the tunnel currents through the left and right-hand 
junctions, Si^i^{uj) and Sij^Ij^^lu) respectively, and the charge noise of the SET island Sq{uj). Using the expression 
for the full current in terms of tunnel currents through the two junctions and the charge conservation relation, we can 
write the current noise spectrum as a sum of the three termstS^ 



Siiiuj) = aSi^i^{Lu) + hSi^i^{uj) - ahuPSqiui), 



where 



— OO 
OO 



dTCOs{ujT) [{QL{t + T)QL{t)) " {QL{t)f 
dTCOs(cJT) [(Qfl(< + T)Qfl(t)} - {QR{t)f 

dr cos{ujT){QD{t + T)QD{t))- 



(11) 

(12) 

(13) 
(14) 



In the last equation, the SET island charge noise, Sq{uj), is identified with the noise in Qnit) (this connection 
will be derived in Sec. V). Notice that the charge noise acts to reduce the full noise spectrum, compared to the 
individual tunnel current noise spectra, reflecting the well-known fact that current noise is suppressed in double 
barrier structures i^S 

The zero-frequency noise of the tunnel currents through the two junctions are equal, because of charge conservation, 
and hence Sij{0) = Sjj^i^{0) = S'/^/^(0). However, in order to obtain the full current noise at finite frequency we 
must calculate each of its three components. The noise in the tunnel currents Si^Ij^{uj) and S'/^/^(ci;) are obtained 
first, using an equation of motion method, then the charge noise, Sq{lu) is calculated using a similar approach. The 
three elements of the current noise are then combined to obtain the full spectrum. 



IV. TUNNEL CURRENT NOISE SPECTRA 



The noise in the tunnel current across the SET junctions can be calculated using an approach due to MacDonalc&Si 
which relies on the fact that the number of charges that have passed across the junction as a function of time can be 
counted. For the SET-resonator system, the description in terms of a probability distribution is easily extended to 
record the number of electrons passing across one or other of the tunnel barriers linking the SET island to the leads. 

Starting with the left-hand junction, we follow MacDonald2S42i in introducing the statistically stationary variable. 



Ql^Ql- {Ql) 



(15) 
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and write, 



where, 



d{Ql{T)) . , ,^ 

OT 



(16) 



QL{t')dt' - {Ql)t 



T-2 2 



(17) 

(18) 
(19) 



with n(r) the number of electrons which passed across the left-hand junction in time r and /q = (/) = (Ql), the 
average current. The MacDonald approach has been used to calculate the current-noise in a number of mesoscopic 
systemsSSiSS and the subtleties of the method are described particularly carefully by Ruskov and Korotkov.^-^ 

The probability distributions for the SET-resonator system are readily generalized to include the number of electrons 
that have passed through cither one of the junctions. In order to calculate the noise in the current through the left- 
hand junction, wc define P^^^^j^j(x, u; i) as the probability that, after time t, n electrons have passed across the 
left-hand junction, there are N{N + 1) electrons on the SET island and the resonator is at position x with velocity 
u. The master equations [Eqs. and 10)] now take the modified form 



dPl 
dt 



2 '^^N 
UJnX- 



du 



fipn 
I— ^ 

dx 



ou ' 



G dP^ 

2m^ du 



dP, 



N+l 



dt 



Lolix - Xo)- 



dP\ 



N+l 



dP\ 



N+l 



du 



dx 



d_ 

du 



leUP, 



N+l 



Tlp;};! - ^rPn 

G dP^^, 



2rii? dt 



i?e2 



N+l 



TrP\ 



R^N 



+ Pn] 
mujQXo 



(20) 



xp^m) 



where P^{x, u; t) = PJ^{x, u; t) + PJ^_^_^{x, u; t). Notice that only the terms which describe the transition from charge 
state iV -I- 1 to iV with an electron tunnelling out through the left junction couple to a change in n. 

The rate of change of the average of the squared number of charges passing through the left-hand junction in time 
T is given byS£ 



dr 



dx I du 



Pn+i{x,u;t) + P^{x,u;t) 



= 2 ri(n) 



JV+l 



muj^XQ 



{xn)i 



N + l -I 



(22) 
(23) 



where we have re-ordered the sum over n in the first term, and substituted from Eqs. (|20|l and (|21|) . The averages 
(. . ■)n+1: arc with respect to the sub-ensemble of systems with iV + 1 charges on the SET island, i.e. 



(...) 



N+l 



E 



dx I du{. . .)PJ^_^_i{x,u;t). 



The initial conditions for the distributions, P^j^^^-^^ (x, u; t — 0), are given by the corresponding steady-state distri- 



butions for the SET-resonator systemji^ pj^}j^^-^Jx,u) 



P. 



^■(Ar+l) — 0) — <5n,0-f^r(jv+l) (-^J ^) 



and hence we have 



for all T. 

The overall expression for the noise in the current through the left-hand junction is given by 



N+l 



K {xn) 



N+l 



Tt 



Xq 



+ elo > sm{ujT)dT. 



(24) 
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The quantities (n(r))Ar+i and (a;n(T))jv+i can be obtained by deriving their equations of motion and those of the 
other average quantities they couple to. The equations of motion for moments of the form (x^u'^n) and (a;^u*n)jv+i 
are given by 

^{xPu'^n) = J2j ^"^J duxPu'^n PJ^^^{x,u;t) + PJ^{x,u;t) , (25) 
{xPu''n}N+i = J2j J duxPu''nP'^+^{x,u;T) (26) 



d 



respectively. Substituting for PJ^^-^^^j^^{x,u;t) from the appropriate master equation [Eqns. 1)201) and H21|l ]. we can 
derive the closed set of equations, 

^{n)r,+i = ^-^TnT-{rL+Tn){n)N+i + ^^^{xn) (27) 
dr e He'' 



— (xn) = (un) + I l{x)n+i 



— {xn) = {un)+rL{x)N+i ^t^{x')n+i (28) 



— (wn) = -uj„[{xn) - xo{n}N+i\ (ux)n+i + Tl{u)n+i - le{un) (29) 

-^{xn)N+i = {un)N+i+Tji{xn) - {TL + Tii){xn)i^+i + ^^^^^^^^{x^n) (30) 

^{un)N+i = [(a;n)Ar+i - a;o(n)Ar+i] + T/jXiin) + ^^^^^(a;un} - (Fl + r;j)(wn)Ar+i - 7e(iin)Ar+i (31) 

^{x'n) = 2{xun) + Tl{x')j,+, - !!^{^^)^^, (32) 

■^{v^n) = -2ujl[{xun) - XQ{un)N+i\ + Tl{u'^)n+i - "^^"f^" (u^a;)Af+i - 2-fe{u'^n) + ^^^t (33) 

-^{xun) = -tJo[(a;^n} - xo(na;)7v+i] + (w^ri) + rL(xu)Ar+i (x^u)Af+i - je{xun). (34) 

The averages which do not involve n are constants (i.e. independent of r) and are evaluated for the steady-state 
distributions, pj^jj^_^^^{x,u)m^ The details of the evaluation of these averages are given in Appendix A. 

The set of coupled differential equations [Eqs. (|2 7134(1 ] is readily integrated numerically to give {n)N+i and {xn)N^i, 
from which the left-hand junction current noise spectrum can be obtained. Since we are interested in the number of 
electrons passing through the left-hand junction in a time t, all quantities involving n will clearly be set to zero at 

T = 0.28 

Integrating the equations of motion for {n)N^i and {xn)N+i and performing the Fourier integral numericallj^ 
(incorporating an appropriate long-time cut-off), we can obtain the current-noise spectrum, Si^i^{uj) for various 
choices of the system parameters. The current noise with the coupling to the resonator set to zero, S^^j^ [uj), can also 
be calculated analytically for comparison,™ 



2/o 



(35) 



(Fl + TrY + cj^ 

Initially, we consider the simplest case, where there is no external damping and the resonator substrate is at zero 
temperature (i.e. 7ert — and = 0), and examine the dependence of the tunnel current noise spectrum on the 
resonator frequency and the strength of the coupling to the SET. Figure |21 shows Sij^i^{llj) for several choices of the 
dimensionless coupling and frequency parameters, k and e.^^ For ease of comparison, the frequency axis is scaled by 
the renormalized resonator frequency, Q = LOTtjd ■ The spectra are dominated by sharp, asymmetric, peaks arising at 

= 0, the renormalized resonator frequency = 1 and at twice that frequency, 17 = 2. None of these three features 
appear in the current noise spectrum for a SET without a resonator. At high frequencies, beyond the range shown in 
Fig. m all of the spectra reduce to the usual single-junction value 2e/o. 

The effect of the resonator on the noise spectrum depends on both the strength of the coupling between the 
resonator position and the SET junction tunnel rates, measured by k, as well as the underlying dynamics of the 
resonator, which in the absence of extrinsic temperature or damping, is conveniently parameterized by 7^X4 = Ace^. 
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FIG. 2: Current noise spectrum of the left-hand junction, Si^i^{ijj), for a range of resonator parameters as a function of the 
scaled frequency = ujTt/e ■ The full spectrum is shown in (a) while the insets (b) and (c) show the details of the peaks at 
n = 1 and 2, respectively. All calculations are performed for 7eTt = = and = Tn. 




O. 



FIG. 3: Comparison of the contributions to the left-hand junction current noise due to the {xn) n+\ and (n) jv+i terms in Eq. 
(I23I I. given by F{u}) and G{ijj) respectively. The curves have been displaced vertically for clarity. From top to bottom, the 
curves are: Si^^i^ (lo) / {2eIo) , F{u!) and G{uj). In each case k — 0.1, e = 0.3, 7eTt = = and F^ = F^j. 



The coupling strength, k, plays a dual role as it controls both the amount by which the resonator motion alters the 
electron tunnel rates, and the resonator dynamics through the intrinsic damping 7^. Fig. [5] shows that the heights of 
the peaks depend sensitively on both k and 7^. The peak heights increase with an increase in k, at fixed 7^, whilst 
an increase in 7, at fixed k (i.e. an increase in e) reduces the heights of the peaks. Notice, however, that the peaks at 
r2 = and 2 reduce more strongly as % increases than the one aX — 1. 

We can understand the origin of the features in the current noise spectrum by splitting it into its two constituent 
parts: that arising from {71)^+1 and that from {xnj^^i. Figure |31 shows a noise spectrum for a particular choice of 
parameters (k = 0.1, e = 0.3), and the two functions F{uj) and G{lo), which add to give the noise spectrum 



2eIo 



F{uj) + G{uj). 



(37) 



The functions F{uj) and G(w) are essentially the Fourier integrals of {xn)N+i and {n)N+i, respectively, with the 
appropriate linear and constant terms subtracted off in each case. Written out explicitly, they have the form 



F(.) = ^ Tdr 
-'0 Jo 



K K 



(xn) 



N+l 



Xq 



sin(Li;T) 



-fo Jo 



dT 



(l-«) 



Io + Tl (e(n) 



AT -1-1 



ar) 



sin(tt>T), 



(38) 
(39) 
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FIG. 4: Effects of finite extrinsic damping and temperature on the current noise through the left-hand junction, Si^i^ (u)), for 
the coupled system with k = 0.05, e = 0.2. Panel (a) shows shows the noise spectra for 7e = (upper curve), 7e = 7i (middle 
curve) and 7e = 57i (lower curve), with = in each case. Panel (b) shows the noise spectra for Q = 0.2 (upper curve), 
B = 0.1 (middle curve) and = (lower curve) with 7^ = 67^ in each case. In both plots the corresponding noise spectrum 
for the uncoupled system is shown as a dashed curve. All calculations are performed for Vl ~ Pjj. 



with 

e (Tl - K/Tt) 

The function G{u!), which contains the (n)jv_|_i term, has a single feature, at the renormahzed resonator frequency. 
In contrast, the function F{uj), which contains the {xn)N+i term, has three features: one at = 0, one at the 
renormahzed frequency, 17 = 1, and one at = 2. We can understand these features in terms of non-trivial correlations 
between the variables x and n. From the spectrum G{uj), we see that the n term has a component oscillating at 
frequency e' and since x will also oscillate at the same frequency, it is not surprising that the combination of the two 
dynamical variables in the term (a;n)jv+i gives rise to a component oscillating at 2e' as well as at e'. The explicit 
dependence on x in {xn)N+i also makes the features in F(lu) more sensitive to the magnitude of 7^ than those in 
G{uj). Hence in the full spectrum the peaks that are due to {xn)N^i alone (at = 0, 2) reduce more rapidly with 
increasing than the one with contributions from {xn)j^^i and {n)i\[^i (at SI = 1). 

The enhancement of the low frequency noise is to be expected because of the effective correlations between successive 
electrons induced by the interactions with the resonator. It is known from other contexts^ that zero-frequency noise 
can be increased when an electron tunnelling event is able to change the internal state of a mesoscopic system (such 
as a quantum dot). In this case, an electron passing through the SET alters the state of the resonator, which in turn 
alters the tunnelling rates for subsequent electrons. However, any excitation of the resonator is damped out over time 
and hence the zero-frequency peak in the spectrum follows a very similar pattern to the other peaks: it increases with 
the coupling strength, k, but reduces with increasing damping, 7^. 

It is the approximation of the gate capacitance by the linear form, Cg{x) ~ Cg(l — x/d), which means that there are 
no further peaks in the current noise beyond that at = 2. If higher order terms had been included in the expression 
for the gate capacitance then they would have lead to higher order terms in the position dependence of the tunnel 
rates, which we would expect to generate peaks at higher multiples of the resonator frequency. However, within the 
regime we consider here, where x/d 1, the additional peaks at higher frequency would be much smaller than those 
at 17 = 1 and 2. 

The effects of extrinsic damping and finite background temperature on the tunnel current noise spectrum are 
illustrated in Fig. 01 with the appropriate noise spectrum for the uncoupled system shown for comparison. It is clear 
that a non-zero value of 7e rounds down each of the peaks at 17 = 0, 1 and 2. Furthermore, the peaks at 17 = 1 and 2 
are also shifted to slightly lower frequencies, reflecting the fact that the extrinsic damping causes an extra reduction 
in the resonator frequency (beyond that due to the interaction with the electrons). The influence of extrinsic damping 
gradually fades at higher frequencies where the spectrum is hardly affected by the resonator anyway. A non-zero 
background temperature counteracts the damping, increasing the height and broadening the peaks. 

The comparison with the uncoupled case in Fig.^ also reveals a small suppression of the noise just after the peak 
at 17 = 1, and a similar feature after the peak at 17 = 2. We can understand the presence of this suppression in a loose 
way by considering the actual motion of the resonator as a competition between two components. Firstly, we can 
consider a resonator coupled to a thermal bath, in this case the mechanical motion would fluctuate strongly leading 
to an enhancement of the noise at the mechanical frequency. "'^^ Of course these thermal-like fluctuations arise already 
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from the interaction between the resonator and the electronsji^ and are increased by the presence of an extrinsic 
bath. However, it seems hkely that the interaction between the electrons and the resonator must also induce some 
coordination between the electronic and mechanical motion which would tend to make the current flow more regular 
and hence suppress the noise, at least at the resonator frequency. As we would expect from this simple picture. Fig. 
^ shows that increasing the background temperature rapidly washes out the regions of noise suppression. 

The noise in the current through the right-hand junction can be obtained in exactly the same manner as that 
through the left-hand junction. Using MacDonald's formula, we have 



K (xn) 



N 



elo > sm{uJT)dT, 



(41) 



where in this case n(r) is the number of electrons to pass through the right-hand junction in time r. 

The master equations [Eqs. Q and are also readily modified to include the number of electrons tunnelling 
(from right to left) across the right-hand junction. A closed set of equations of motion that allow (n)jv and {xn)^ 
to be determined by numerical integration are obtained following the same procedure as that used for the left-hand 
junction. 

Carrying out the numerical integration and then the Fourier integral, as before, we can compare the resulting noise 
spectrum through the right-hand junction with that through the left-hand junction. As would be expected for the 
very symmetrical SET layout we consider here, the current noise spectra of the two junctions are very similar, but 
away from $1 = they are not identical. In fact, the magnitudes of the peaks in the left and right-hand tunnel noise 
spectra at fi = 1,2 differ by up to a few percent for the parameters used here. The differences in the two spectra 
arise because the resonator breaks the overall symmetry of the circuit. 



V. CHARGE NOISE ON THE SET ISLAND 



Apart from the contributions due to the tunnel currents across the two SET junctions, the full current-noise 
spectrum also contains a contribution arising from fluctuations in the charge stored on the SET island, 

/oo 
dTcos{ujT){QD{t + T)QD{t)). (42) 
'OO 

The average charge on the SET island is a constant (i.e. (Qd) = 0), thus we can use MacDonald's formula to show 
explicitly that this term is indeed the charge noise spectrum, 

Sq{lu) = 2a>-i /"^^%^sin(c^r)dr (43) 
Jo 

poo 

= 4/ dTCOs{u;T)[{QDit + T)QD{t))-{Ql)], (44) 
Jo 

where an integration by parts has been performed. 

Within our model of the SET, the excess charge can only be either zero or one electron, hence — —e{QD)- 

Since [Qoit)) is simply the average occupancy of the SET island multiplied by the electronic charge, — e, we can use 
the master equations to write down a closed set of coupled equations that link the average excess charge, position, 
{x), and velocity, (u), of the resonator, 

J^{QD{t)) = -eTn ^{Qnit)) - "^{xit)) (45) 
j^{xit)} = {nit)) (46) 

j^iuit)) = -u;',{xit)) - '^{Qoit)) ~ iMt)). (47) 

Note that these equations are equivalent to the mean-coordinate equations derived in Ref. 

Since the equations of motion for the probability distribution are Markovian, the equations of motion for two-time 
correlation functions of the charge can be obtained from the equations for the average quantities [Eqs. H45l47|l ] using 
the regression theoremt^ Thus we obtain, 

^{Quit + T)QD{t)) = -eYn{QD(t)) - ^{QD^t -f T)Qu{t)) - ^^^^{x{t + r)g,,(i)) (48) 
dr en en 
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Q 

FIG. 5: Frequency scaled charge noise spectrum, uj'^Sq{lo), for a range of resonator parameters, with jeTt = = and 
Fl = Ffl. 

■^{xit + T)QDit)) = {uit + T)QD{t)) (49) 

^{u{t + T)QD{t)) = -iolixit + T)QD{t)) - "^{Qoit + r)QD{t)) ~ le{u{t + T)QD{t)). (50) 
dr e 

The initial conditions (at r = 0) are given by the steady state properties of the system. The charge noise spectrum is 
then obtained by integrating the coupled equations to obtain {Qoit + T)QD{t)), and then performing the necessary 
Fourier integral. 

The frequency scaled spectrum, uj'^Sq{ll)), is shown in Fig.[Slfor a range of resonator couplings and frequencies. The 
charge noise spectrum of the SET island, like the function G{uj) considered earlier, contains only a single peak, at 
the renormalized resonator frequency. This is not surprising as the equations of motion from which Sq{lu) is obtained 
[Eqs. H48l50|l ] have a structure which is very similar to those which determine G{llj) [Eqs. H27l29|l ]. Figure 13 shows 
that the height of the fl = 1 peak in (jj'^Sq{u!) is almost independent of k and e (provided k and e are much less than 
unity, and also non-zero), this is because the peak height in the charge noise itself is approximately determined by 
the ratio of the coupling to the intrinsic damping K/(7iTt) = 



VI. FULL CURRENT NOISE SPECTRUM 



Once the noise spectra of the tunnel currents and the charge on the island have been calculated, the full current- 
noise spectrum is readily obtained by combining its components in the appropriate way [as given by Eq. (|11|) ]. Figure 
shows the full current-noise spectrum for the particular choice of the parameters a — b. The overall spectrum is 
dominated by peaks ai Q — 0,1 and 2, and largely flat elsewhere. The peak in the charge noise spectrum suppresses 
the peak in the tunnel noise spectra at = 1, so that the peak at = 2 can become larger in the overall spectrum. 
For a combination of sufficiently strong coupling and weak damping, the zero-frequency noise exceeds the Poissonian 
value 2e/o. 

In the limit where k goes to zero, the full current noise spectrum is readily obtained analytically. In this limit the 
left and right-hand junction current noise spectra are the same, and the full current noise (spectrum assuming a = b) 
is 



1 



1 - 



Y2 



r2 



(51) 



where F = F^, + F/?. This is of course the usual result for sequential tunnelling through a symmetric double-barrier 
structure. ^^■^■^ Thus in strong contrast to the coupled system, the current noise of the uncoupled, symmetric (Tb. = Fl), 
SET would be flat: S'J^(tj)/(2e/o) = 0.5. 

Figure |H1 shows in detail the dependence of the spectrum on the resonator coupling, k and intrinsic damping 
7iTt = Ke^. Once the frequency of the spectra is scaled by the renormalized frequency, a relatively simple behavior 
emerges. Increasing the coupling strength increases the height and widths of the peaks, whilst increasing 7^ (for fixed 
coupling) simply rounds the peaks down without affecting their widths. Although the heights of the peaks at f2 = 1 
and 2 scale with coupling in a very similar way in the individual tunnel current spectra (see Fig. [21 , the relative 
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FIG. 6: Full current noise spectrum, Sii{u}), for a range of resonator parameters, assuming a — b with = 0, 7eTi = 0. The full 
spectrum is shown in (a) while the insets (b) and (c) show the details of the peaks at f2 = 1 and 2, respectively. All calculations 
are performed for Fl = Tr. 




(a) 



I 



0,54 




FIG. 7: Full current noise spectrum including extrinsic effects. The influence of extrinsic damping at 9 = are shown in (a) 
where the curves (from top to bottom) are for 7^ = 0, 7i and 67^. The effect of a finite background temperature is shown in (b) 
where the curves (from top to bottom) are for = 0.2,0.1 and 0, and 7^ = 67^ in each case. The calculations are performed 
for K = 0.05, e = 0.2, a = 6 and Fl = Fr. 



magnitudes of the peaks at f2 = 1 and Q = 2 in the full spectrum are very sensitive to the magnitude of n. This is 
because the peak in the frequency-scaled charge noise uj'^Sq{uj) is almost completely insensitive to changes in k and 
hence gives rise to a proportionately stronger suppression of the f2 = 1 peak for smaller k values. 

For the relatively low frequencies shown in Fig. the noise is higher than the uncoupled value, even away from 
the immediate vicinity of the peaks. This increase in the background noise level depends on k, but not e, and occurs 
because the resonator breaks the symmetry of the tunnel rates through the two junctions. Only for the uncoupled 
system do the physical tunnel rates correspond to the parameters F^ and F^. When the coupling to the resonator is 
switched on, the tunnel rates for both junctions depend linearly on the position of the resonator, but with opposite 
sign, and the strength of this coupling is given by k. The coupling to the resonator can thus be thought of as giving 
rise to an effective renormalization of the tunnel rates and since the noise is a minimum for the symmetric tunnel 
rates of the uncoupled system, F^ = F^ (which we choose here), the coupling always gives rise to an increased 
background noise level. Although it is not shown in the figures, the current noise reduces to the uncoupled value at 
large frequencies, f2 ^ 1, since any asymmetry in the junction tunnel rates becomes unimportant at large frequencies. 

The influence of extrinsic effects on the current noise spectrum are illustrated in Fig.[7| As with the tunnel current 
spectra, extrinsic damping reduces the peak heights and the background noise level, whilst increasing the background 
temperature has the opposite effect. Furthermore, it is interesting to note that extrinsic effects can alter the relative 
heights of the peaks at = 1 and 2. 

The full noise spectrum could equally well be calculated for a variety of cases where a / 6, reflecting the balance of 
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capacitances in the system more accurately, as discussed in Sec. II. However, now that we have obtained the individual 
spectra of the noise in the tunnel currents and the charge noise, it is clear that for small changes of the coefficients 
about the values a = b = 1/2, there would be no qualitative change in the overall current spectrum and so we have 
not investigated such cases explicitly. 

VII. CONCLUSIONS 

The presence of a voltage-gate consisting of a nanomechanical resonator has a dramatic effect on the current-noise 
spectrum of the SET. The low frequency noise of the SET can be substantially enhanced and two additional peaks 
appear in the noise spectrum: at the frequency of the resonator (renormalized by the interaction with the SET) and 
at twice that frequency. Without extrinsic damping and with the background temperature set to zero, the peaks 
in the current noise are sharp, with heights controlled by the the strength of the resonator-SET coupling and the 
magnitude of the intrinsic damping constant. 

The features in the current-noise spectrum reflect the fact that the coupled dynamics of the SET and the resonator 
lead to important correlations between charges flowing through the SET and the position of the resonator. In 
particular, the interaction between the resonator and the electrons on the SET island provide a mechanism by which 
the passage of one electron through the system can affect the motion of subsequent electrons. Hence it is not surprising 
that the noise in the measured current does not conform to any simple picture based on a combination of the usual 
noise spectrum of the SET and the mechanical noise spectrum of a damped harmonic oscillator. 

The calculation presented here is for the particular case of a classical resonator coupled to a SET as a voltage gate. 
Very recently a number of other authors have considered aspects of the current noise in a similar nanoelectromechanical 
system known as the electron shuttle^ in both the classical^^*^ and quantum regimes. ^4 In particular, the current 
noise spectrum of a classical shuttle system obtained by Isacsson and Nord^ shows peaks at the frequency of the 
mechanical system and its first harmonic in close correspondence with the spectrum of the SET-resonator system 
described here. 
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APPENDIX A: CALCULATION OF STEADY-STATE PROPERTIES 

The moments of the steady-state distribution which appear in the set of equations of motion describing the motion 
of charges through the SET [Eqs. 127I34|| ] can be obtained analytically from the master equations [Eqs. and Q]- 
Adding Eqs. and Q leads to an equation of motion for u; t) which in turn can be used to obtain an equation 
of motion for the moment A similar equation of motion for the moment {x'^u^)i^^i can be obtained directly 

from Eq. |2Jl . 

Assuming a steady state, the equations of motion for {x'^u^) and {x"'u'^) lead to the following linear equations, 

2m^ 

k^ij^„ [(x"+iu^-i)Ar+i - xo{x^u^-')m+i] = - (fc7e + + Th)(u''x'')n+i (A2) 

+k{k - 1)^{u''-^x-)n+i + Tniu'^x") + !!^{x"+^u''), 

respectively. By considering these equations for appropriate values of k and n, a set of linear equations is generated 
which can be solved to obtain the necessary steady-state moments, assuming only that the steady-state probability 
distribution is an even function of the velocity, u. 
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